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Abstract 

CL We investigate invasions from a biological reservoir to an initially empty, heterogeneous habitat in the presence of advection. The 
Q habitat consists of a periodic alternation of favorable and unfavorable patches. In the latter the population dies at fixed rate. In the 
" ^ former it grows either with the logistic or with an AUee effect type dynamics, where the population has to overcome a threshold to 
lP grow. We study the conditions for successful invasions and the speed of the invasion process, which is numerically and analytically 
^"'investigated in several limits. Generically advection enhances the downstream invasion speed but decreases the population size 
1 'rjf the invading species, and can even inhibit the invasion process. Remarkably, however, the rate of population increase, which 
j quantifies the invasion efficiency, is maximized by an optimal advection velocity. In models with AUee effect, differently from the 
logistic case, above a critical unfavorable patch size the population localizes in a favorable patch, being unable to invade the 
habitat. However, we show that advection, when intense enough, may activate the invasion process. 
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Invasions of alien species are widespread phenomena, in 
principle affecting every ecosystem, usually with dramatic 
consequences on the native community, constituting a ma- 
jor threat to biodiversity (Vitousck ct al, 1997; Mooney, 
2000; Pimentel et al, 2000). At the scale of interest for man- 
agement purposes, i.e. the geographic scale, invasive species 
move across a heterogeneous landscape characterized by 
favorable and unfavorable areas. The presence of abiotic 
heterogeneity, in fact, characterizes most of natural habi- 
tats and plays a key role in invasion processes, influencing 
their rate of spread and outcome (Shigesada and Kawasaki, 
1997; Hastings et al, 2005; Melbourne et al., 2007). 

Alongside with the empirical interest for the problem, 
several modeling efforts have been dedicated to the un- 
derstanding and prediction of the spatial spread of invad- 
ing organisms in heterogeneous environments. Within the 
framework of reaction diffusion models, building on the pi- 
oneering theoretical works of Skellam (1951) and Kierstead 
and Slobodkin (1953) on the "critical patch size" problem, 
Shigesada et al (1986) gave a seminal contribution consid- 
ering the invasion (propagation) of a population through a 
periodic heterogeneous environment (see also Weinberger, 



2002; Kinezaki et al, 2003). The problem was extended in- 
cluding advective transport to study persistence and propa- 
gation of passively dispersing populations in oceans (Mann 
and Lazier, 1991; Abraham, 1998) or rivers (Speirs and 
Gurney, 2001; Pachepsky et al, 2005; Lutscher ct al, 2006). 
The importance of the interplay between heterogeneity and 
advection has been recently reviewed by Ryabov and Bla- 
sius (2008). Moreover, the role of both advection and land- 
scape spatial structure is clearly relevant also to the disper- 
sal of plants (Hastings et al, 2005), whose seeds are trans- 
ported by winds. 

In this article, we focus on the interplay between abi- 
otic heterogeneity and advection in invasions. We describe 
the dynamics in terms of an advection-reaction-diffusion 
model, which allows for mathematical tractability and 
quantitative predictions, e.g., on the spreading rates. 

We consider an infinite system where a population sta- 
bly saturates the carrying capacity on one side of the sys- 
tem and possibly invades the remaining part of the envi- 
ronment, which is assumed to be heterogeneous. Our set- 
ting is quite general and widely applicable. In particular, it 
is relevant to situations in which one has a practically infi- 
nite biological reservoir of a species invading an empty ter- 
ritory characterized by abiotic heterogeneity. For instance, 
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the above setting may be relevant to situations in which 
invasions can suddenly become possible for the removal of 
a climatic barrier due to climate changes (Mooney, 2000). 
Another relevant case is when a species stably populating 
a lake invades an effluent characterized by a certain degree 
of heterogeneity and stream velocity. This is one of the key 
early-stage processes related to the spatial control of inva- 
sions in lakes' networks (Havel et al, 2002). Other examples 
concern the spreading of wind-pollinated plants in a hetero- 
geneous environment (Davis et al, 2004) or spores carried 
by the wind (Kot et al, 1996). 

More specifically, the habitat consists of a periodic alter- 
nation of unfavorable and favorable patches, as in Shigc- 
sada et al (1986). The population dies at a fixed rate in 
unfavorable regions, and grows in favorable ones according 
to either a logistic or an Allee effect dynamics. We are in- 
terested in determining the conditions for invasions to be 
possible and in understanding how invasion speed and effi- 
ciency depend on the mechanisms at play. 

With the logistic dynamics, in the absence of advection, 
this problem was pioneered by Shigesada et al (1986), while 
Lutscher et al (2006) considered both advection and hetero- 
geneity in reference to the "drift paradox" problem (Spcirs 
and Gurney, 2001). Going beyond these works, we find 
asymptotic expressions for both the invasion speed and the 
rate of increase of the population size. The latter quantity 
essentially estimates the rate at which the number of invad- 
ing individuals grows and, thus, provides a suitable measure 
of the efficiency of the spreading process. Indeed, especially 
in invasive species control, it is important to quantify the 
potentiality of growth of an alien population, and not only 
the speed at which it colonizes the territory. We anticipate 
that , remarkably, larger invasion speeds do not necessarily 
imply more efficient invasions. 

The logistic case (decreasing per capita growth rate) is 
then contrasted with the case of positive density depen- 
dence corresponding to a demographic Alice effect (Allee, 
1938; Dennis, 1989), which accounts for a reduced repro- 
ductive power at low densities. The importance of Allee ef- 
fects for the invasion and control of non-native species was 
emphasized by Taylor and Hastings (2005) and Tobin et al 
(2011). It is interesting to mention that even in homoge- 
neous habitats the presence of the Allee effect can decrease 
the invasion speed or even halt the population spreading if 
the initially occupied area is too small (Lewis and Kareiva, 
1993) (see also Vercken et al (2011) for recent field obser- 
vations). We find that the interplay between heterogeneity 
and advection becomes very subtle in the presence of the 
Allee effect. In fact it may happen that a persisting popula- 
tion, unable to invade new territory, becomes able to spread 
in the presence of strong advection. This effect should be 
taken as a cautionary note from the standpoint of con- 
trolling invasive species, telling us that advective transport 
should be considered. For instance, after strong weather 
events, or in regions characterized by prevailing winds, ne- 
glecting the effects of advection could lead to the erroneous 
prediction of a population unable to invade, whereas it ac- 



tually propagates over the territory. 

The material is organized as follows. In Sect. 2 we present 
the model, and in Sect. 3 we qualitatively discuss its phe- 
nomenology. Sections 4 and 5 present and discuss the main 
results on invasions with the logistic and the Allee effect 
model, respectively. Finally, in Sect. 6 we summarize the 
results. 

2. Model 

The evolution of the population, 9(x,t), is governed by 
the advection-rcaction-diffusion equation 

d t e + vd x e = Ddle + f(e,x). (i) 

The diffusion coefficient D and the advection velocity v are 
assumed to be constant. Heterogeneity is introduced in the 
growth term, f(9, x), which depends on the position x. The 
habitat consists of a periodic alternation of unfavorable 
and favorable patches of sizes l u and If, respectively. In 
the elementary cell [0 : C] (where L = l u + If denotes the 
spatial period), we take 

f gje) < x < l u 
f(6,x) = { ya( ' ~ . (2) 

( g f (0) i u <x<c 

In the unfavorable regions the population is assumed to die 
at a constant rate r u , so that g u (9) = —r u 9 (r u > 0). In the 
favorable regions we consider two dynamics. The first is the 
classical logistic model (with carrying capacity normalized 
to one) 

g f (6)=r f 9(l-e), (3) 

Tf being the intrinsic growth rate. Equation (1) with the 
logistic term but without advection was firstly studied by 
Shigesada et al (1986). Recently, Lutscher et al (2006) in- 
cluded advection focusing on the "drift paradox" problem. 

Secondly, accounting for a positive correlation between 
population density and per capita growth rate at small den- 
sities — the Allee effect (Allee, 1938; Dennis, 1989) — , we 
consider the threshold model 

g f (e)=r f m^{(e-e c )(l-e),0}, (4) 

prescribing that the population grows only when 8 > 8 C 
(otherwise it stays constant). Notice that (4) recovers 
(3) for 8 C = 0. We remark that the model (4) represents 
an intermediate case between weak and strong Allee ef- 
fect (Courchamp et al, 2008, see also Sect. 5 for further 
discussions). To the best of our knowledge, models with 
Allee effects have been mostly investigated in homoge- 
neous habitats (Petrovskii and Li, 2003). In heterogeneous 
habitats we are aware of only a few studies with integro- 
diffcrence models incorporating different dispersal kernels 
(see, e.g., the recent work by Dewhirst and Lutscher, 2009; 
Pachcpsky and Levine, 2011). 

We now specify the settings in which Eq. (1) is studied. 
We consider model (1) with boundary condition 9(0, t) = 1, 
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Fig. 1. Population evolution with logistic dynamics and boundary condition (BC) 6(0, t) = 1 (a,c) or with periodic BC, 6(0, t) = 6(nL, t) with 
n = 3(b,d). Panel (a): Successful invasion for l u =2, J* = 1.8, e = l and u = 0.7. The solid curves represent 9(x,t) at three successive times 
while the dotted curve corresponds to the case without advection (u = 0), for comparison. The inset shows J°° dx6(x,t) vs t. Panel (b): 
Using the same parameters as in (a) but with periodic BC. The dashed rectangle corresponds to the initial condition, while the gray thick 
line is the population at stationarity. The black line superimposed on the gray one is obtained taking a cut of the interior of the traveling 
front from (a). Panel (c): Unsuccessful invasion obtained in the same settings as in (a) but for If = 1.4, notice that here for u = the 
invasion is still possible, meaning that advection is responsible for the halt of the invasion process. The three solid curves refer to the same 
time instants of (a), while the dotted one is the u = case. Inset as in (a). Panel (d): parameters as in (c) but with periodic BC. Now the 
population goes extinct. The inset displays 9(x,t) at three successive times (from top to bottom) showing that 6 — > exponentially. 



mimicking the case in which on the left of the origin (x < 0) 
the population constantly saturates the carrying capacity, 
while the population is initially absent in the x > region, 
i.e. 0(x, 0) = for x > 0. With this choice for the boundary 
conditions the invasion process must be considered from 
left to right (i.e. from the biological reservoir at x < to the 
positive real axis) . In this case depending on the sign of the 
advection velocity we can consider (downstream) invasions 
with the flow (i.e when v > 0) or (upstream) invasions 
against the flow (i.e. when v < 0). 

It is useful to formulate the model in non-dimensional 
variables. To this aim we exploit known results about the 
logistic growth model without advection, namely for the 
standard FKPP equation (Fisher, 1937; Kolmogorov et al, 
1937). The FKPP equation develops traveling fronts char- 
acterized by the propagation speed v = 2wDrf and width 
£o = \jDjrf. It is then natural to measure lengths in units 
of £o , time in units of the inverse growth rate in the favor- 
able patches 1/r/, and the advection velocity in units of vq. 
We thus define the non-dimensional variables x' — x/£o, 
t' = trt, u = v/vq. Dropping the primes, Eq. (1) made 
non-dimensional reads 

d t + 2ud x = d 2 x + f{6,x). (5) 

The factor 2 in the advection term results from our choice 
to fix u — 1 as the non-dimensional propagation speed in 
the homogeneous FKPP system. We can now introduce 
e = Tu/rj which is the death over growth rate ratio, and 



lf,u = £f,u/£o which are the non-dimensional sizes of the 
patches (L = £/£o = l u + If). In this way, with reference 
to Eq. (2) we have g u (0) = —^0 and 

g f (0) =0(1-0), (6) 

for the logistic model, while with an Allee effect it becomes 

5/ (0)=max{(0-0 c )(l-0),O}. (7) 

3. Model phenomenology 

We now present the basic phenomenology of the model, 
discussing also the main differences between logistic and 
Alice effect growth models. A successful invasion implies the 
development, far from the boundary, of a traveling front, 
characterized by a stationary and spatially periodic bulk 
(Fig. la). In such a case the total population in the in- 
vaded habitat, J Q dx0(x,t), asymptotically increases lin- 
early with time (inset of Fig. la). Faster growing popu- 
lations mean more effective invasions. Conversely, Fig. lc 
shows a typical case of unsuccessful invasion: no traveling 
front develops and the total population remains bounded 
in the limit of long times (compare the insets in Figs, la 
and lc). 

For positive advection velocities, the problem of identi- 
fying the conditions for successful invasions is directly re- 
lated to determining under which conditions Eq. (5) with 
(6) in a finite system (of size nL with n integer) with pe- 
riodic boundary conditions (0(0, t) = 0(nL,i)) admits a 
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Fig. 2. Population evolution with the Allee effect model (7). The black and gray curves are obtained without advection (u = 0) and with 
advcction (u = 0.5), respectively. The other parameters are e = 1, It = In = 8 and 9 C = 0.001. The different panels refer to: (a) BC 0(0, t) = 1, 
(b) the same of (a) but in logarithmic scale, (c) periodic BC with n = 3 elementary cells. The dashed line in (c) is the initial condition. 



non-vanishing stationary solution, starting from a generic 
non-zero initial condition. This is clearly shown by the per- 
fect superposition of the bulk of the traveling front with 
the stationary solution of the periodic boundary condi- 
tion problem (Fig. lb). The reason for this link is that the 
bulk region of the traveling front (which is stationary and 
spatially periodic) satisfies the same boundary value prob- 
lem of the finite system with periodic boundary conditions 
(BC). Therefore, to determine whether downstream inva- 
sions are successful it is enough to study whether persis- 
tence is possible in the finite system with periodic BC. Con- 
versely, an unsuccessful invasion implies that in the finite 
system the population goes extinct, exponentially in time 
(inset of Fig. Id). 

As we will see in the following, for small patch sizes, the 
qualitative behavior of the different growth models is very 
similar: invasions benefit from larger favorable regions and 
their speed increases accordingly; moreover, the presence 
of advection enhances the downstream invasion speed but 
decreases the population size of the invading species, even- 
tually halting the invasion (see Fig.l). Conversely, for large 
enough unfavorable patches, dramatic differences appear. 
In the Allee effect model, the population can persist in the 
absence of advection but localized in a region of finite size, 
if it initially occupied that area (see Fig. 2c), being un- 
able to invade new territories. This is quite different from 
the logistic model where invasions and persistence are al- 
ways linked. Even more striking is the role of advection. 
Figures 2a, b show that suitable values of the advection ve- 
locity can activate the invasion of an otherwise localized 
population. 

In the following we present the results for the logistic 
and Allee effect model separately as the level of analyti- 
cal understanding is quite different. The possibility to use 
the linear analysis framework in the logistic model, indeed, 
allows us to systematically derive the conditions for inva- 
sions and asymptotic expressions for the invasion speed and 



efficiency. This approach cannot be used for the Allee ef- 
fect model, which is studied mainly numerically and with 
heuristic arguments. 



4. Results for the logistic growth model 

4.1. Persistence in a closed periodic system 

In this section wc focus on the conditions for persistence 
in a closed system with periodic BC, which correspond to 
those for successful downstream invasion. Moreover, when 
the population is able to persist, we study how its size be- 
haves as a function of the advection velocity. It is worth 
noticing that the closed system setting is interesting also in 
consideration of recent experiments where bacterial popu- 
lations are grown in heterogeneous conditions (Dahmen et 
al, 2000; Lin et al, 2004; Perry, 2005). 



4.1.1. Critical patch size and critical advection 

Starting from a population different from zero in a single 
favorable patch, with periodic BC, Eq. (5) with (6) admits 
cither the trivial solution 6 — 0, meaning that the habitat 
is unable to sustain the population, or an asymptotically 
stationary non- vanishing solution, when persistence is pos- 
sible. To determine the conditions for the latter, it is suffi- 
cient to identify when the solution 9 = becomes linearly 
unstable, as briefly sketched in A. 

For l u and e fixed it is possible to show the existence 
of both a critical size of the favorable patch, If, such that 
extinction occurs if If < If , and a critical advection velocity 
u c such that for a range of values of If the system goes 
extinct if u > u c . The implicit relation between the critical 
values reads 
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Fig. 3. Extinction (white) and survival (gray) regions in the plane 
(lf,u). The boundary curve between the two regions gives the de- 
pendence of the critical velocity u c on If . The two vertical lines pin- 
point the values l*„ and el u . The Roman numbers on top label the 
three regions discussed in the text. Data refer to l u = 2 and e = l. 
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The above equation was also found by Lutscher et al (2006) 
using boundary conditions different from ours. Notice that 
for \u\ > 1 the term yl — u 2 ' becomes imaginary, and in 
Eq. (8) the identities sin(iz) = isinh(z) and tan(iz) = 
i tanh(z) must be employed. We also remark that Eq. (8) 
is left unchanged by the substitution u — > — u due to the 
symmetries of model (5) with periodic BC. Therefore, we 
can limit the analysis to u > 0. 

Concerning the critical size of the favorable patch, in the 
case of small unfavorable regions (l u -C 1) the advection 
has not a great effect: at the leading order, ps el u , i.e. 
the dependence on u is negligible. However, for large un- 
favorable patch sizes, we have (see also Speirs and Gurney 
(2001); Ryabov and Blasius (2008)) 
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showing that advection worsens the survival conditions. 
This effect can be deduced noticing that the advection term 
changes the growth/death rate into 



e(x) — > e(x) 



(9) 



where e(x) is the spatially dependent growth rate, tak- 
ing the values 1 and — e in the favorable and unfavorable 
patches, respectively: essentially the effective growth rate 
is decreased by u 2 while the effective death rate is increased 
by the same amount. Equation (9) can be derived from (5) 
via the transformation 8(x,t) — > 8(x, t)e ux (Dahmen et al, 
2000; Ryabov and Blasius, 2008). However, this transfor- 
mation changes the value of the density at the boundaries 
making the solution of the periodic BC case more cumber- 
some. 



When If > If, the population can be driven to extinction 
by intense advection, exceeding a critical velocity u c that 
can be computed from Eq. (8). Figure 3 shows in gray the 
region in parameter space (lf,u) where the population is 
able to survive, for e and l u fixed. We can identify three 
regions (as labeled on the top of Fig. 3): (I) for If < If 
the population goes extinct for any value of u; (II) for < 
If < survival is possible below a critical velocity u c ; 
for If — > el u we have that u c — > oo; (III) for If > el u the 
average growth rate is positive and the population survives 
for any value of u. 

4.1.2. Effects of advection velocity on the population size 

Now we study how the size of the population depends on 
the advection velocity in order to characterize the transition 
from survival to extinction and to derive some results to 
be used later (Sect. 4.2). In particular, we are interested in 
the behavior of the average biomass defined as 



B = lim(0(x,i)) = lim i / 0(x,t)dx , 

t^OO t— ¥OQ ±j J 



(10) 



in the limit of large times, when the solution is stationary. 
In the above expression, thanks to the periodicity of 9{x, t), 
we considered the biomass present in an elementary cell. 
Given the habitat properties, the biomass is a function of 
the advection velocity u, B = B{u). 

With periodic boundary conditions, B is an even function 
of u, B(u) = B(—u), so that, assuming that B(u) is a 
smooth function, for u -C 1 we expect 



B(u) = B(0) [1 



nil 



0(u 4 )] , 



(11) 



with r) some positive constant, as confirmed by numeri- 
cal simulations (dotted curves in Fig. 4). The biomass de- 
creases with u because the net effect of advection is to in- 
crease/decrease the death/growth rate, as from Eq. (9). 
Equation (11) agrees with results of Dahmen et al (2000) 
for the linearized dynamics. Even though they also claim 
that with the complete equation non analytic behaviors 
(i.e. B(0) — B{u) oc |u|) may appear due to the nonlinear- 
ity However, our simulations always confirmed (11). 

The behavior (11) holds both in region II and III of Fig. 3. 
In region III, where survival is possible for any value of the 
advection velocity, we have that in the limit u — > oo the 
average biomass attains a finite limiting value given by (see 
B) B(u) = A/l f + 0{A/u 2 ) , where A = lf-el u , is nothing 
but the average growth rate times L. Curve (a) in Fig. 4 
shows B{u) at the transition between region II and III, i.e. 
for A = 0. In this case u c = oo, so that B(u — > oo) = and 
the function 



B{u) = 



B(0) 
1 + Pu 2 



(12) 



provides a very good fit of B(u) for any value of u. Notice 
that Eq. (12) implies (11) with n = f3. 




Fig. 4. Biomass B vs u with l u = 2 and e = 1. Curve (a) corresponds 
to lj = el u = 2, at the border between region II and III; the other 
curves refer to values of If in region II: If = 1.8 (b) and If = 1.5 (c). 
Symbols show the results of simulations while the solid curves shows 
the fitting functions: for (a) Eq. (12) with = 0.726; for (b) and 
(c) Eq. (13) with u c = 1.823, /3 = 0.626 and u c = 0.794, /3 = 0.494, 
respectively. The theoretical values for u c obtained from (8) are 1% 
close to the fitted values. The dotted parabola shows the behavior 
(11). Inset: B(u) vs u c — u for (b) and (c) close to the transition. The 
straight lines show the slopes 2B(0)/(u c (l + obtained from 

(13) when u — > u c . 

In region II, the critical velocity u c is finite and, by defi- 
nition, B(u) — > when u — > u c . As typical in phase transi- 
tions, we should expect B{u) ~ (u c — u) u for u c — u <C 1, 
with v some exponent characterizing the extinction transi- 
tion. Assuming a smooth behavior it is reasonable to expect 
v = 1, as confirmed by the inset of Fig. 4 and supported by 
analytical approaches by Dahmen et al (2000) valid in the 
limit l u — > oo. Finally, assuming the simplest functional 
form consistent with the symmetries and regularity prop- 
erties of B(u) we end up with the expression 



B(u) = B(0) 



1 - (u/u c y 



(13) 



which is consistent with (11) for u«l giving n = f3 + u~ 2 , 
and with (12) for u c — > oo. Curves (b) and (c) in Fig. 4 show 
B(u) for two values of If within region II. For both values 
one can observe the very good agreement between numeri- 
cal data and Eq. (13). The above results show to what ex- 
tent advection decreases the average population size, and 
provide a characterization of the advection-induccd popu- 
lation extinction. 

4.2. Effects of heterogeneity and advection on the invasion 
speed and efficiency 

As discussed in Sect. 3, Eq. (8) also provides the condition 
for downstream invasions (i.e. when u > 0) from a reservoir 
(on the left) to a heterogeneous habitat (on the right, as 
in Fig. la). As for upstream invasions (i.e. when u < 0), it 
is necessary to understand when Eq. (5) admits solutions 
which develop a periodic traveling front advancing with a 
positive speed u p , for long times and far from the boundary. 



In the following, we show how the speed u p can be derived, 
and discuss the general conditions for invasions. 

Assuming that, far from the boundaries, a traveling front 
develops, following Shigcsada et al (1986) we can write 
0(x, t) = Q(z)g(x), where z — x — 2u p t accounts for prop- 
agation with velocity u p (the factor 2 deriving from our 
choice of the non-dimensional variables, see Sect. 2). The 
function Q(z) describes the traveling front modulated by a 
periodic function g(x) = g(x + L) due to the habitat peri- 
odicity. For the computation of the invasion speed we can 
use the linearized dynamics assuming that the traveling 
component has an exponential leading edge &(z) cx e~ sz . 
As detailed in C, we end up with an implicit relation be- 
tween the invasion speed u p and the shape parameter s of 
the traveling front (see also Lutscher et al (2006)) 



T>(u p , s;u) — cosh(go_L) — cosh(q u L u ) cosh(g/Z/) 

r,2 



'In 



2q u qf 



smh(q u l u ) smh(q f lf) 



(14) 



where go 



\J u 1 + e + 2su p and 



<1J 



\Jv? — 1 + 2su p . From (14) one derives u p (s; u) and the in- 
vasion speed can be obtained by computing min s {u p (s; u)}. 
As far as we know, there is no analytical expression for the 
minimum, and numerical computations must be employed. 
In the sequel, with some abuse of notation we will denote 
with s(u) the value of s for which the minimum is realized 
and with u p (u) the minimal speed. 
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Fig. 5. Regions where invasions are successful (gray) and unsuccessful 
(white) in the plane (lf,u). For u > we have downstream invasions. 
For u < 0, upstream invasions are possible only for — 1 < u < 0, see 
text for explanation. Data refer to the case l u = 2 and e = 1. 

The gray area in Fig. 5 displays the region in the plane 
(lf,u) where invasions are possible. Such a region was nu- 
merically determined by solving Eq. (14) and finding the 
values of // and u for which the propagation speed u p (u) = 
min s {u p (s; u)} exists and is positive. For u > the curve, 
separating white and gray regions, coincides with that de- 
rived from Eq. (8) (shown in Fig. 3). For u < 0, it ap- 
proaches the asymptote u = — 1 for large If. Indeed for 
very large favorable patches the system should recover the 



homogeneous habitat result u p = 1 + u, so that invasions 
are impossible for u < —1 (Lutscher et al, 2006). 

We now focus on the invasion speed u p for downstream 
invasions {u > 0). The upstream case was considered in 
details by Lutscher et al (2006). Simulations (not shown) 
suggest that for the invasion speed behaves linearly 

in u, i.e. 



u p (u) ~ u p (0) + au . 



(15) 



both in region II and III. The above result holds also for 
small negative u. To determine a the first step is to expand 
T>(u p , s;u) in powers of u. At the first order the expansion 
yields 



V(u p , s; u) = V{u v , s; 0) + f(s)u + 0{u 2 ) = , 



(16) 



with f(s) — Lsinh(Ls). For small u the minimum of (16) 
is realized at u p (u) = u p (0) + 5u p and s(u) = s(0) + Ss, 
with Su p ,Ss ~ 0(u). Expanding now Eq. (16) in Su p and 
Ss, one gets 



SUr, 



M0)) 



5 Up D(u p , S ;0)| Up(0)i 



«(0) 



au . 



(17) 



From the numerical values of u p (0) and s(0), we computed 
the value of a obtaining a perfect agreement with simula- 
tions. 

We now consider the behavior of the invasion speed for 
large advection velocities. When I f is chosen in region II the 
invasion is halted for u > u c . In region III, the population 
can invade the habitat for any u and the invasion speed 
approaches another linear behavior for u 3> 1. In Fig. 6, 
we contrast two cases: when If = el u and If > el u , inside 
region III. 




Fig. 6. Invasion speed u p vs u for If = 2 (black curve) and If = 2.5 
(gray curve) , together with the linear behavior u p = u (dotted line) . 
Inset (a): u p — u vs u for If = 2 (black curve) and If = 2.5 (gray 
curve). Inset (b): u p — u vs If (symbols) compared with the asymp- 
totic prediction (18), approximated in the simulation using u = 100. 
Data have been obtained with l u = 2 and e = 1. 

As discussed in Sect. 4.1.2, when If = el u , the average 
biomass vanishes (B(u) —> 0) for u u c = oo. It is thus 



reasonable to expect that the contribution to the invasion 
speed comes only from advection. Therefore, asymptoti- 
cally we expect u p — > u as shown in Fig. 6, though the 
convergence of u p — u to zero can be rather slow (Fig. 6a) . 
Conversely, inside region III, u p ~u reaches a finite value for 
u — >• oo (Fig. 6a). As heuristically derived in D, for u ^> 1 
one expects 



i p {u) -u = A/If , 



(18) 



where A — If — el u . Strictly speaking, the above result 
holds in the limit of u — > oo and If — > oo but, as shown in 
Fig. 6(b), it is in fairly good agreement with the numerical 
results also for finite values of If. 

To summarize, in region III the invasion speed u p is well 
approximated by two different linear behaviors 



Up[u) ~ u p (0) + au for small u 
Up{u) ~ A/If -\-u for large u 
where a is given by Eq. (17). 



(19) 



4.2.1. Invasion speed in rapidly and slowly varying 
environments 

We now study the invasion speed in two limiting cases 
for which some analytical results can be obtained, namely 
when the habitat is finely fragmented (l u j <C 1) or subdi- 
vided in large patches (l u j S> 1 with 7 = lf/l u fixed). 

When lf u -C 1, expanding Eq. (14) at the lowest order, 
one obtains the explicit expression 



u p (s; u) 



1 

2~s 



(20) 



Retaining higher order terms, it is possible to show that the 
dependence of u p on u is linear up to the fourth order in 
lf, u , where a term proportional to u 2 appears. Minimizing 
Eq. (20) yields 



i p {u) = y / A/L' + 



(21) 



as also obtained with homogenization techniques (Lutscher 
et al, 2006). Equation (21) shows that for finely fragmented 
habitats the intrinsic propagation speed is as in the homo- 
geneous habitat once the growth rate is substituted with 
the average growth rate A/L, as found by Shigesada et al 
(1986) for u = 0. Moreover, Eq. (21) implies that for l u j — > 
0, a in Eq. (15) is equal to unity. We finally observe that 
Eq. (21) holds only for A > 0, that is in region III of Fig. 3. 
Indeed for l u j — > region II shrinks to zero. 

In the limit of very large patch sizes, lf_ u — > 00, from 
Eq. (14) we obtain 



q L = q u l u + q f lf . 



(22) 



where q — it+s, q u = y /, u 2 + e+2su p and 9/ = \/u 2 — l + 2su p ' 
(see C). It is interesting to compute the limit maintaining 
the ratio 7 = lf/l u constant so that Eq. (22) becomes 
<7o(l +7) = q u + 79/, which has the explicit form 



1 




Fig. 7. Invasion speed u p (0) vs l u with 7 = lf/l u fixed and e = 1. 
Continuous lines correspond to the numerical solution of Eq. (14) for 
7 = 1 (black) and 7 = 1.5 (gray). Dashed lines correspond to the limit 
values for rapidly and slowly varying environments, respectively. The 
former is zero for 7 = 1. Symbols show u p measured in numerical 
simulations with the Allee effect model (7) with 6 C = 0.001 (see 
Sect. 5). Black and gray symbols correspond to 7 = 1 and 1.5, 
respectively. The vertical bar marks the critical unfavorable patch 
size (27) above which invasions are impossible. 



(l+7)(w + s)=y / w 2 +e+2u p s' +7 a/" 2 - l+2u p 



(23) 



The above formula does not depend on l u and / f separately 
but only on their ratio 7, meaning that for l u j — > 00 the 
propagation speed u p approaches a limit value that depends 
only on 7 and e. The limits 7> 1 and 7 <C 1 are quite trivial 
and consistent with intuition. In both cases, neglecting sub- 
leading terms, squaring both sides of (23) one finds u p (s). In 
the former limit, If 3> l u (negligible unfavorable patches), 
the homogeneous result u p (u) = 1 + u is retrieved. In the 
latter, l u ^> If (negligible favorable patches), the condition 
for the minimum of u p (s; u) is realized for imaginary values 
of s, meaning that invasions are not possible. 

In the special case 7=1 and u = 0, Eq. (23) becomes 

from which it is easy to 



2s 



derive the invasion speed: 



x/2 



u p s + e + 



1 



u p (0) 



l + e 2 + (l-e)\/l + e + 



1 



2VT 



3/2 



(24) 



Notice that, once the proper correspondence between nota- 
tions is made, the result (24) coincides with that obtained 
by Hamel et al (2010) using a different technique. Further 
specializing to the case of equal growth and death rates 
(e = 1) Eq. (24) reduces to u p (0) = 2 1 / 2 3~ 3 / 4 . 

In the presence of advection with u<l, the computation 
for 7 = 1 and e = 1 can be easily extended obtaining 
u p (u) = 2 1 / 2 3~ 3 / 4 + (2/3)u that is Eq. (15) with a = 2/3. 
For u> 1, we were unable to obtain analytical results, but 
we expect the phenomenology discussed in Fig. 6 to apply. 

In Fig. 7 we show the invasion speed u p , obtained by 
numerically solving (14) for u = 0, at varying l u j with 7 = 
lf/lu = 1 and 1.5. We also show the asymptotic values for 
luj — > and l u j — ¥ co. In rapidly varying environments 



u p (0) = y/ (j — e)/(l + 7)' while in slowly varying ones 
u p (0) goes to the finite value (23) which for 7 = 1 is given 
by (24). 



4.2.2. Efficiency of the invasion process as a function of 
the advection velocity 

The invasion speed u p measures the velocity of popula- 
tion advancement. We now focus on the rate of increase 
of the population size (inset of Fig. la), which provides a 
measure of the efficiency of the invasion process. The suit- 
able quantity to look at is the rate of increase of the total 
biomass 



B r {u) 



j t J 6(x,t)dx, 



which is the slope of the curve shown in the inset of Fig. la. 




Fie 



Rate of increase of the biomass B r (u) vs u for e = 1, l u = 2 



and If = 1.5 (boxes) and If = 1.8 (circles) (see also Fig. 4). 

Figure 8 shows that there is an optimal advection velocity 
which maximizes the rate B r (u). So that, in spite of the 
fact that larger advection velocities imply larger invasion 
speeds, the efficiency of the invasion process is maximized 
at a specific value u* of the advection velocity. This means 
that even though for u > u* the invasion speed increases, 
the number of invading individuals decreases, which implies 
a less effective invasion. 

We now provide a heuristic argument to explain the ori- 
gin of an optimal advection velocity. At stationarity, the 
population, advancing at constant speed u p (u), increases 
its size at a rate 



B r (u) ~ B(u)u p (u) 



(25) 



where B(u) is the average biomass (10). For u <C 1, B{u) is 
well described by Eq. (11), i.e. B(u) decreases quadratically 
with u. On the other hand, as from Eq. (15), the invasion 
speed increases linearly with u. Using the above considera- 
tions and Eq. (25), we obtain that the increase in invasion 
speed will dominate at very small u while the quadratic de- 
crease of B(u) will dominate at larger u, producing the bell 
shaped behavior observed in Fig. 8. 



The above argument is based on a low order Taylor ex- 
pansion which, in principle, could cease to be valid for val- 
ues of u at which the maximum of B r (u) is attained. We 
numerically found that the behavior reported in Fig. 8 is 
general and that, typically, the maximum of B r (u) is real- 
ized for values of u for which the Taylor expansion is still 
a valid approximation. 

4.3. Discussions 

Shigesada et al (1986) have shown that even when the 
average growth rate is negative (i.e. A = lf — el u < 0) there 
exists a critical size of the favorable patches above which 
a population can invade new territories. We have shown 
that the main effect of advection is to increase the critical 
size l*f for the invasion to be possible. In particular, there 
always exists a critical advection velocity u c above which 
no invasion is possible, unless the average growth rate is 
positive. 

As for the invasion speed, we recover the results obtained 
for u = by Shigesada et al (1986), and for u ^ by 
Lutscher et al (2006), who focused on upstream propaga- 
tion (i.e. in our setting u < 0). However, here, we mainly 
focused on the downstream invasion speed, i.e. when the 
advection velocity favors the invasion (u > 0). Our analy- 
sis shows that, provided u < u c , advection always increases 
the invasion speed. In particular, both for small and large 
advection velocities u the invasion speed is linear in u but 
with different prefactors (see Eq. (19)). In some interesting 
environmental limits, moreover, we analytically computed 
the invasion speed. Our results agree with those found with 
different techniques by Hamel et al (2010), and are inter- 
esting in view of the ensuing discussion on Alice effects. 

Our results on the dependence of the population size (the 
biomass) on advection extend similar ones derived by Dah- 
men et al (2000) in the limit of large unfavorable patches. 
Although the invasion speed is enhanced by advection, the 
biomass decreases at increasing the advection velocity. As 
a consequence, a faster invasion speed does not necessarily 
imply a more efficient invasion process. Indeed, the suit- 
able quantity to judge about the effectiveness of the inva- 
sion process is the rate of increase of the biomass. We found 
the remarkable new result of an optimal advection veloc- 
ity maximizing such rate and, hence, the efficiency of the 
invasion process. This maximum originates from the oppo- 
site role of advection on the invasion speed u p and on the 
biomass B(u): the former increases linearly with u while 
the latter decreases quadratically with u. The balance be- 
tween these two behaviors leads to a maximum in the in- 
vasion efficiency. 

5. Results for the Allee effect growth model 

Some features of the logistic model seem quite unreason- 
able from an ecological point of view. For instance, as shown 
in Fig. 7, the invasion speed approaches an asymptotic value 



when the sizes of the patches are enlarged holding fixed 
their ratio, regardless of the size of hostile regions (i.e. even 
for /„ — > oo). Intuition would suggest that extremely large 
unfavorable regions should slowdown, and eventually sup- 
press, invasions, as observed, e.g., in soil organisms (Bai- 
ley et al, 2000). A slowdown of the invasion process should 
be expected, indeed, anytime there is a positive correlation 
between population density and per capita growth rate, i.e. 
in the presence of Allee effects (Allee, 1938; Dennis, 1989). 
In general, one speaks of Allee effect when for small densi- 
ties the growth rate is negative — strong Allee effect — or 
positive, but smaller than for larger density — weak Alice 
effect (Wang and Kot, 2001). Relevant works and their re- 
lations with our problem are discussed in the next section 
(in particular, sec Table 1). 

We consider the Allee effect threshold model (7), that we 
recall here 

g f {6) = max{(0-0 c )(l-0),O}. (26) 

This model corresponds to a situation in between the strong 
and weak Allee effect and is convenient because it reduces 
to the logistic model (6) for 6 c — 0, easing the comparison. 
At the end of this section we will briefly discuss different 
Alice effect models. 

Studying the problem without advection allows us to 
identify the main consequence of the Alice effect, namely 
the existence of a critical unfavorable patch size /* above 
which invasion is impossible for any size of the favorable 
region and the population can persist localized as shown 
in Fig. 2. At stationarity, in the bulk of the traveling front 
we have that, if 9 m is the value at the beginning of an un- 
favorable patch, then at the end of the unfavorable region 
the density will reach the value 9 m = 9m cxp (— \fel u ). 
Given the growth term (26), for the population to prop- 
agate we must require that 9 m > 9 C (see also Dcwhirst 
and Lutscher (2009) for a similar argument applied to an 
intcgro-differcncc model), which implies the inequality 

''^Tf^ilt) (27) 

for the unfavorable patch size. In general, 9m cannot be 
estimated analytically; however, setting 9m = 1 gives a rea- 
sonable upper bound. 

The existence of /* is evident from Fig. 7 where symbols 
denote the results of numerical simulations obtained with 
the growth model (26), holding constant 7 = and 
increasing the patch sizes. For small sizes the qualitative 
behavior of the model (26) is similar to that of the logis- 
tic model: u p increases with L. A dramatic difference ap- 
pears at large sizes: for the logistic model the invasion speed 
reaches an asymptotic value, while for the Allee effect one 
it decreases and, eventually, the invasion process is halted 
when l u ~ lui regardless the size of the favorable patch. 

In Figure 9, we show the behavior of the system without 
advection in the plane (l u ,lf), for two values of 9 C . As al- 
ready discussed, for the logistic growth model there exists a 
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Fig. 9. Regions of the plane (l u ,lf) characterized by localization, 
propagation and extinction of the population with the Allee effect 
model (26) with u = 0, e = 1, C = 0.001 (black) and 8 C = 0.01 
(gray). For each value of 9 C , two curves separate the three different 
possibilities as shown in the inset with shaded areas for the case 
Be = 0.001. The vertical lines in the main figure show /* obtained 
from (27) with 8m = 1. The dotted curve displays the dependence 
of I* for the logistic model (6). 

critical favorable patch size If above which the population 
can survive and invade new territories. The critical value If 
remains finite for l u — » oo as found by Ludwig et al (1979); 
Shigcsada et al (1986) and also in this paper. With the 
Allee effect growth model (26) the phenomenology is differ- 
ent and more interesting. When the size of the unfavorable 
patches is small (i.e. l u -C 1) the system essentially behaves 
as the logistic model: the unfavorable patch is so small that 
typically the decrease in density will not cause the popu- 
lation to fall below 9 C . From a quantitative point of view, 
l*f is slightly larger than the logistic value (the effect being 
more pronounced for larger 9 C ). though this cannot be fully 
appreciated from Fig. 9 due to the scale. The main qual- 
itative change with respect to the logistic case manifests 
when the size of the unfavorable patches approaches the 
value i* given in Eq. (27). For unfavorable patches at least 
this large, propagation becomes impossible for any size of 
the favorable patches, even though the population does not 
necessarily go extinct. Indeed, as highlighted in the inset of 
Fig. 9, a new region in the plane (/«,?/) appears, where the 
population can persist locally but cannot propagate (see 
Fig. 2): it localizes in a single favorable patch (if initially it 
was in that patch) . For this localization regime to exist it is 
necessary that I / is wide enough to sustain the population, 
as theoretically derived in the homogeneous case by Lewis 
and Kareiva (1993) and found in field data by Vercken et 
al (2011). 

Therefore, without advection but with the Allee effect 
growth term (26), we have that the population goes extinct 
if If < If, propagates if If > If but l u < I*, and localizes in 
a single favorable patch when this is large enough to sustain 
the population (If > If) but the unfavorable patch is too 
large to allow the propagation, i.e. /„ > £*. 

We now discuss the effects of advection. For small un- 




Fig. 10. Invasion speed u p vs u for the Allee effect model (26) with 
8 C = 0.001, e = 1 and l u = lf = 8, i.e. within the localization region 
of Fig. 9. Gray symbols correspond to values of u for which the 
population remains localized, while black ones to values for which 
it propagates. Inset (a): rate of increase of the biomass B r vs u 
obtained using the same parameters of the main figure. Inset (b): 
invasion speed vs If for u = 0.2 and l u = 8. At increasing If the 
population passes from extinction (empty symbols) to localization 
(gray symbols) and finally to propagation (black symbols). The last 
transition is possible only thanks to advection. 

favorable patch sizes, i„ < 1, since the model with Alice 
effect behaves quite similarly to the logistic model, also 
the effect of advection on the dynamics is very similar be- 
tween the two models. A critical advection velocity u c ex- 
ists, above which the population goes extinct. However, for 
larger unfavorable patches, while in the absence of advec- 
tion propagation is inhibited and the population remains 
localized (light gray region in the inset of Fig. 9), sustained 
advection induces the remarkable qualitative changes ob- 
served in Fig. 2. In Fig. 10 we show a numerical measure- 
ment of the invasion speed u p as a function of u for a habi- 
tat with unfavorable and favorable patch sizes, l u and // 
respectively, chosen in the localization region of the system 
without advection. As it can be seen, while the population 
remains localized at small advection velocities, it suddenly 
becomes able to propagate invading the whole environment 
when the velocity of the medium becomes large enough. An 
intuitive explanation for such a behavior is that, thanks to 
advection, the population can now travel through the un- 
favorable patch more rapidly, finally reaching the next fa- 
vorable patch with a density above the threshold 9 C , i.e., 
advection enhances the value above which no propa- 
gation is possible. In fact, at stationarity, in the unfavor- 
able region it is easy to see that the density behaves as 
9(x) = 9m exp[(u — vtlr Therefore, the same argu- 
ment which lead to Eq. (27) now yields 

g(u)= * log(^). 
ye + u A — u \ t>c / 

The above formula predicts that l^(u) grows with u, so 
that even if in a medium at rest the population is localized, 
i.e. /„ > lu(0), in advective media there will be a value 
of the advection velocity such that l u < /*(m), allowing 
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Fig. 11. Advection-induced invasions for different Allee effect models with C = O.OO1, l u = S, If = 7 and e = l. (a) With u = the population 
evolving with the Allee effect models as labeled is unable to invade, (b) With u = 0.5 invasion is possible but for the cubic model (28). As for 
model (29) we used b= —1. (c) Comparison of the per capita growth rates gf(9)/0 in the three considered models together with the logistic 
case gf(8)/8 = 1—9 (gray line). Notice that for comparison purposes all models have been normalized in a such a way that maxg{gf(8)/8} = 1. 



the population to propagate, as shown in Fig. 10. Once 
invasion is permitted by advection, the propagation speed 
and the rate of increase of biomass behave similarly to the 
same quantities in the logistic model (compare Fig. 10 main 
figure with Fig. 6 main figure, and Fig. 10a with Fig. 8). 

The results shown in the inset (b) of Fig. 10 further illus- 
trate the importance of advection in such model. There, we 
show the measured propagation velocity u p as a function of 
If, for constant advection u and keeping l u > I* (0) fixed, so 
that in the absence of advection the population would be 
localized even if I f > lj . As one can see two transitions are 
observed: from extinction to localized, mainly due to the 
increase of If; from localized to propagating, which is pos- 
sible only thanks to advection. It is interesting here to note 
that something similar happens if we consider a spatially 
dependent diffusion coefficient, as done by Shigesada et al 
(1986) for the logistic growth model. In particular, with- 
out advection (u = 0) and denoting with 8 = D u /Df the 
ratio between diffusivities in the unfavorable and favorable 
patches, Eq. (27) modifies in I* w ^Jdje \i\{9 M / 9 c ) ■ This 
result tells us that a larger diffusivity in the unfavorable re- 
gions (5 > 1) allows the population to propagate also with 
larger unfavorable patch sizes. Therefore, if we fix the size 
of the unfavorable patch, either increasing the diffusivity or 
the advection, the residence time in the unfavorable patch 
decreases leading to a smaller depletion in population den- 
sity. 

Through this section we limited the numerical analysis to 
the Allee effect model (26) which is intermediate between 
the case of weak and strong Allee effect. It is thus worth 
to test the robustness of the above findings by considering 
different Alice effect models. In Figure 11 we show the evo- 
lution of invasions obtained with three different models of 
the Alice effect, cither in the absence (Fig. 11a) or in the 
presence (Fig. lib) of advection when l u > I*. In particu- 
lar, we compare model (26) with two models of Alice effect, 



namely the standard cubic term used to model strong Alice 
effects (Lewis and Kareiva, 1993; Pctrovskii and Li, 2003; 
Almeida et al, 2006) 



g f (9) = 8(1 - 8)(9 - 9 C ) , 



(28) 



and a slight modification of the model introduced by Wang 
and Kot (2001) 



9f(9) 



be 



i 



< 9 < 9 C 

e c < 9 < i 



(29) 



The latter can model either strong (b < 0) or weak (if < 
b < 1 and 1 — 9 C > b) Alice effect. The marginal case b = 0, 
for small values of 9 C , is practically indistinguishable from 
the Alice effect model (26). Figure 11a shows that with u = 
invasions are unsuccessful in all the considered models. 
For model (29) we only show the case b = — 1. However, the 
result (not shown) holds also for positive (not too large) 
values of b, meaning that also with weak Allee effects the 
population can be localized. Advection (Fig. lib) is able 
to activate the invasion process in all cases but for the cu- 
bic model (28). The origin of such a difference is readily 
understandable from Fig. 11c where we show the density- 
dependent per capita growth rate gf(9)/9 in the above 
models compared with the logistic case (the gray line) . As 
one can see, the growth rate for the cubic model, although 
positive, remains very small for a large interval above the 
threshold value 9 C . As a consequence, even if the population 
reaches the favorable region with 9 > 9 C , it cannot propa- 
gate. In the other cases the growth rate is large enough to 
allow the population to grow and propagate. 

Concluding, the phenomenology of the advection- 
induced invasion is rather general in the presence of (either 
weak or strong) Allee effect, but critically depends on the 
growth rate realized close to the threshold value 9 C . In 
perspective, it would be interesting to investigate what 
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invasion, homogeneous diffusion+advection Lewis and Kareiva (1993), 
persistence, (short+long range) Petrovskii and Li (2003), 

localization Almeida et al. (2006) 



persistence heterogeneous diffusion 

(short range) 



Shi and Shivaji (2006) 



invasion, heterogeneous diffusion+advection This work 
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invasion, heterogeneous dispersal kernel Dewhirst and 

localization, (long range) Lutscher (2009), 
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Table 1 

Models with Allee effects related to the present work. The models in the upper part of the table are advection-reaction-diffusion models; those 
in the lower part of the table are integro-difference models. Further details about models of Alice effects in invasion dynamics can be found 
in tables 3.5 and 3.6 in Courchamp et al (2008). Our work aims at investigating conditions for persistence and invasion in heterogeneous 
environments with advection. 



habitat/advection conditions ensure successful invasions 
for different kinds of Alice effects. 

5.1. Discussions 

To put in perspective our results on the Allee effect model 
it is useful first to briefly recall some known results from 
the literature. The study of Allee effects has received con- 
siderable interest in the past, with growing modeling efforts 
in recent years, when its importance for invasions and con- 
servation issues have started to be well recognized (Taylor 
and Hastings, 2005; Courchamp et al, 2008; Tobin et al, 
2011). For convenience of the reader, we summarize rele- 
vant works related to our problem in Table 1. 

In the case of homogeneous environments and without 
advection, most studies have shown that the success of 
an invasion depends not only on the initial density, but 
also on the size of the initially occupied area (Lewis and 
Kareiva, 1993; Lewis and van den Driessche, 1993; Kot et 
al, 1996; Wang and Kot, 2001). Asymptotic rates of spread 
are typically reduced (see, e.g. Lewis and Kareiva, 1993), 
mostly because of the interplay between dispersal mech- 
anisms and the reduced reproductive power at low densi- 
ties. In discrete-space models (e.g., in Kcitt et al, 2001) 
the interesting phenomenon of population localization (also 



called range pinning) occurs. The presence of advection 
can lead to nontrivial results when combined with density- 
dependent migration (Petrovskii and Li, 2003; Almeida et 
al, 2006). 

In the case of heterogeneous habitats, with few excep- 
tions (see, e.g., Shi and Shivaji (2006) who extended the 
critical patch size problem to the case of weak Allee effects) , 
most works focused on integro-difference models (Dewhirst 
and Lutscher, 2009), sometimes accounting also for the 
effect of discreteness of the population (Pachcpsky and 
Levine, 2011). These studies agree on the fact that the 
presence of Allee effects combined with habitat fragmen- 
tation generally penalizes the success of invasions or, at 
least, slows them down. In particular, crucial for the suc- 
cess of the invasion process is the size of bad patches (De- 
whirst and Lutscher, 2009) with respect to the dispersal 
range. Moreover, the results of these models depend also on 
the choice of the dispersal kernel, which is typically poorly 
known (Hastings et al, 2005). 

Our results show that logistic and Alice effect models 
qualitatively display the same features both in the pres- 
ence and in the absence of advection when the unfavor- 
able patches are not too large, including the nontrivial ex- 
istence of an optimal advection velocity maximizing the in- 
vasion efficiency. However, in the Allee effect case, unlike 



the logistic model, if the unfavorable patches become too 
large invasions can be halted. When the Alice effect is in- 
cluded, in fact, if both l u and If are large, the population 
can persist locally, but is unable to invade other patches, 
since it cannot cross unfavorable regions without being too 
severely damped. This means that while in the logistic case 
the thresholds for persistence and for invasions coincide, 
these are in general different when Allee effects are present. 
The same observation was made in the context of integro- 
diffcrence models (Dewhirst and Lutscher, 2009), and it is 
substantiated also by field observations (Bailey et al (2000), 
see also Hastings et al (2005) and references therein). 

Advcction, however, can alter this picture: if its inten- 
sity exceeds a threshold value, the invasion process can be 
activated, provided that a large enough growth rate is real- 
ized at the population density when entering the favorable 
patches. This peculiar effect of advection, which as far as 
we know was not previously put into light, has obvious eco- 
logical implications for invasive species management strate- 
gies. For instance, the idea to induce Allee effects to control 
the invasion of alien species, which can be implemented in 
several ways (Tobin et al, 2011), should be pursued with 
extreme care in the presence of advection. 

In summary, our results on persistence and invasion in 
the presence of Allee effects extend previous works done in 
the framework of reaction-diffusion models (with homoge- 
neous habitat) and positively correlate also with the results 
of intcgro-difference models, emphasizing the subtle role of 
advection which was missed in that kind of models. 

6. Conclusions 

In this paper we focused on the role of advection on in- 
vasions in heterogeneous environments characterized by fa- 
vorable and unfavorable patches. 

On the one hand, we have shown that in the presence 
of habitat heterogeneity sufficiently intense advection can 
halt the invasion process. Moreover, we argued that the 
efficiency of the invasive process is properly quantified in 
terms of the rate of increase of the invading population. In 
particular, we found that the latter is maximal at interme- 
diate values of the advection velocity. 

On the other hand we have shown that in the presence 
of Alice effects, advection may be beneficial to the invasion 
process turning a persistent but non invading population 
into an invading one. 

An important aspect in evaluating biological invasions, 
which has not been considered in our work, is related to the 
discrete nature of a population, which is made of individu- 
als (Durrett and Levin, 1994; Okubo and Levin, 2001). We 
expect that considering also discrete effects may add fur- 
ther nontrivial effects due to demographic stochasticity, in 
particular close to the transition between successful and un- 
successful invasions. Works in these direction have started 
to appear (Snyder, 2003; Pachcpsky and Levine, 2011). It 
would be interesting to study the effect of advection also in 



the presence of demographic stochasticity. 
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Appendix A. Conditions for persistence in a 
periodic system 

Due to the habitat periodicity, and with the periodic BC, 
we can limit the analysis to the unit cell [0 : L). Denoting 
with d a perturbation around the solution 9 = 0, we have 
that "d is ruled by the linearized version of Eq. (5) 

d t >d + 2ud x ti = dl§ + e(x)§, (A.l) 

where e(x) = —e in the unfavorable patches and e(x) = 
1 in the favorable ones. In the linear analysis framework, 
at leading order, one expects that i3 = e At ip(x) so that 
population extinction is a stable solution if A < and 
an unstable one if A > 0. When unstable, the asymptotic 
solution will be a stationary and spatially periodic solution 
ip{x) as in Fig. la, otherwise if A < it decays exponentially 
as -&(x,t) = -0(a;)e~' A !* (see inset of Fig. lc). Plugging i9 = 
e At ip(x) into (A.l) we obtain the characteristic equation 
for the stationary state 

d^tp - 2ud x ip + (e(x) - A)ip = , (A.2) 
whose general solution is given by 

i> u (x) = A?e xX ^ + A%e* x ™ 0<x<l u 
ip f (x) = A{e xxi " + A f 2 e xX * f) l u <x<L 1 

where 

A^ 2 } =u± Ve + A + u 2 ' = u ± 6„(A) 
A^ 2 } =u± i\J\ — A — it 2 ' = u± ibf(K) 

are the eigenvalues associated to Eq. (A.2). For a solution 
to exist it is sufficient to impose the continuity of densi- 
ties ip u j(x) and fluxes J u j(x) = [2u — dx^ujix) at the 
boundaries between favorable and unfavorable regions. No- 
tice that the conditions on fluxes are required when using 
space-dependent diffusion coefficient, as, for example, in 
Shigesada et al (1986) and Lutscher et al (2006). For con- 
stant diffusion coefficient, as here, it is enough to impose 
the continuity of derivatives of ip. 

Using Eq. (A. 3) and imposing the aforementioned conti- 
nuity conditions we obtain a linear system for the four con- 
stants A^'l. Requiring that this linear system has a non- 
trivial solution we obtain 



cosh(uL) — cos(6/(A)Z/) cosh(b u (A)l u ) 



bl(A)-b}(A) 
26 U (A)6/(A) 



i(b f (A)l f )smh(b u (A)l u ) 



(A.5) 



The largest value of A solving the above equation deter- 
mines the stability properties of the solution 9 = 0. In par- 
ticular, fixing the values of l u , e and u, it is possible to show 
that A > whenever If > I*., where l* f solves Eq. (A.5) 



with A = (Shigesada et al, 1986; Nagylaki, 1975), i.e. 
Eq. (8). For If < 1% the population always goes extinct. 
An equivalent critical value exists for the advection veloc- 
ity, see main text. Let us notice that with u = 0, the above 
condition reduces to that found by Shigesada ct al (1986), 
i.e. -y/F tanh (yfel u /2) = tan(Z//2). For l u <C 1, the above 
equation tells us that ZJ w el u . Notice that the average 
growth rate is given by (If — d u ) / L and that for If > el u 
survival of the population is guaranteed by the fact that 
the average growth rate is positive. As the size of the un- 
favorable patches grows the critical (favorable) patch size 
approaches the limit value l*f = 2 arctan yft corresponding 
to the result of Ludwig ct al (1979). The case e — > oo, i.e. in- 
finite mortality, corresponds to the KISS critical patch size 
I) = 7T (Skellam, 1951; Kierstead and Slobodkin, 1953). 



Appendix B. Derivation of the asymptotic 
expression for biomass. 



Appendix C. Derivation of the dispersion relations 

For long times, as shown in Fig. 1, the bulk of the trav- 
eling front is a periodic function, 9(x, t) = 9(x + L,t + T) , 
where L coincides with the habitat spatial period and T 
is the temporal period; the invasion (propagation) speed is 
then given by L/T. With the chosen BC the propagation 
proceeds in the positive x direction. To derive the prop- 
agation speed we can write 6(x,t) = Q(z)g(x) with z = 
x—2u p t = x+L — 2u p (t+T) and g(x) = g(x+L) being a pe- 
riodic function which modulates the traveling front (Shige- 
sada et al, 1986). We can now take Q(z) cx e~ sz meaning 
that, apart from the periodic modulation g(x), the leading 
edge is exponentially decaying. Plugging the above expres- 
sions in the linearized Eq. (5) yields the equation for g 



g"-2(u + s)g'+[e(x) + 2(u-u p )s + s 2 ]g = 0. 



(C.l) 



If e(x) = 1, g is constant and the homogeneous case is re- 
covered. The linear equation (C.l) has the general solution 



9u(x) 



and gf(x) 



(/), 



.co- 



in the unfavorable and favorable patches, respectively and 



jJL{l = u + s± Ju 2 + e + 2su p = qa ± q u 



= u + s±\/u 2 -l + 2su p = q ± qf . 



At stationarity, for «> 1 we can disregard the diffusive 
term in Eq. (5) obtaining: 



-e0/(2u) < x < l u 

9(l-9)/(2u) l u <x<L 



(B.l) 



where the prime represents the derivative with respect to 
x. Then, denoting with 9m = 0(0) and 9 m = 9(l u ) the 
values of 9 at the beginning of the unfavorable and favorable 
regions, respectively (which correspond to the maximum 
and minimum realized values of 9), Eq. (B.l) is solved by 



9(x) 



^exp(-g) 



, exp 



' X-l-g 
v 2u 



l-fl m [l-exp(^)] 



< x < l u 



l u < x < L 



(B.2) 



Now imposing the periodicity 9m = 9(L) and noticing that 
Qm = 6m exp(— el u /(2u)) we find that 



exp(^)-l 



ex P (t) ~ ex P 1 



(B.3) 



where A = If — el u . The above expression provides a mean- 
ingful solution only for A > 0, i.e. in region III where ex- 
tinction never takes place. Moreover, for u — > oo Eq. (B.3) 
can be expanded to show that it reaches a finite limit 9m = 
A/Z/(l + 0(l/u)). Integrating (B.2) in [0 : L\ one obtains 
an explicit expression for B(u), which expanded for large 
u gives B(u) = A/If + 0(A/u 2 ) . 



Imposing the continuity conditions as in A, we find four 



equations for the coefficients 



and requiring the exis- 



tence of a nontrivial solution, we find that the dispersion 
relations (14) must be satisfied. Notice that (14) was also 
found by Lutscher et al (2006) using different BC. More- 
over for u — Eq. (14) reduces to the equation found by 
Shigesada ct al (1986). 



Appendix D. 
advection 



Invasion speed in the limit of large 



The values of u p (u) — u and 9m (the maximum value of 
the population density) can be written in general as u p (u) — 
u = f(u, e, If, l u )u and 9 M = g(u, e, If, l u )9 , where u = 
1 and 9q = 1 are the values in a homogeneous environment. 
Taking the ratio between the above equations we obtain: 

f(u,e,l f ,l u ) 

U p (u) -U = — : 7—-9 M ■ 

g(u,e,lf,l u ) 

Now, if we conjecture that the decrease of the difference 
u p (u) — u with u is mainly controlled by the dependence 
of 9m on u, we can make the strong assumption that 
f(u,e,lf,l u )/g(u,e,lf,l u ) w 1, which is true at least for 
If — >• oo as both functions tend to 1. Using Eq. (B.3) of B 
in the limit u — ¥ oo, which implies 9m = A/Z/, we finally 
obtain u p (u) — u = A/ If . This equation is true in the limit 
of u — > oo and If — > oo but also for finite values of If it is 
in fairly good agreement with the numerical results (see 
Fig. 6b). 
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